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The first stars in the Universe are predicted to have been much more massive 



than the Sun 1 3 . Gravitational condensation accompanied by cooling of the 
primordial gas due to molecular hydrogen, yields a minimum fragmentation scale 

O ■ 

of a few hundred solar masses. Numerical simulations indicate that once a gas 

CO ' 

clump acquires this mass, it undergoes a slow, quasi-hydrostatic contraction 
without further fragmentation 1,2 . Here we show that as soon as the primordial 
gas — left over from the Big Bang — is enriched by supernovae to a carbon or 

> 

oxygen abundance as small as ~ 0.01—0.1% of that found in the Sun, cooling by 
singly-ionized carbon or neutral oxygen can lead to the formation of low-mass 

X. . . . . . 

stars. This mechanism naturally accommodates the discovery of solar mass 
stars with unusually low (10 -5 3 of the solar value) iron abundance but with a 
high (10 13 solar) carbon abundance. The minimum stellar mass at early epochs 
is partially regulated by the temperature of the cosmic microwave background. 
The derived critical abundances can be used to identify those metal-poor stars in 
our Milky Way galaxy with elemental patterns imprinted by the first supernovae. 

The microphysics of cooling due to molecular hydrogen (H2) within the first gas clouds 
to condense in the Universe, introduces a characteristic temperature (T c ~ 100-200 K) and 
hydrogen density (n c ~ 10 4 cm -3 ) into the formation process of the first stars 1 ' 2 . The 
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temperature floor is set by the energy separation of the lowest-lying rotational levels of H 2 , 
and the characteristic density leads to the thermalization of these levels, at which point 
cooling becomes less efficient. The minimum mass scale for fragmentation is approximately 
given by the Jeans mass, Mj ~ 4OOM (T/2OO K) 3 / 2 (n/10 4 cm~ 3 )' 1 / 2 . The final mass of a star 
is determined by the accretion process onto the nascent protostellar core, resulting in M* ~ 
aMj with a <^ 0.5 5 ' 6 . In contrast to the formation mode of massive stars (population III) at 
high redshifts, fragmentation is observed to favor stars below a solar mass (population I and 
II) in the present-day Universe. The transition between these fundamental modes is expected 
to be mainly driven by the progressive enrichment of the cosmic gas with heavy elements 
(or 'metals'), which enable the gas to cool to lower temperatures. The concept of a 'critical 
metallicity', Z crit , has been used 7-11 to characterize the transition between population III 
and population II formation modes (where Z denotes the mass fraction contributed by all 
heavy elements). Previous studies 7,8 ' 11 have only constrained this important parameter to 
within a few orders of magnitude, Z crit ~ 1CT 6 — 10~ 3 Z Q , under the implicit assumption 
of solar relative abundances of metals. This assumption is likely to be violated by the 
metal yields of the first supernovae (SNe) at high-redshifts, for which strong deviations from 
solar abundance ratios are predicted 12-17 . The cooling rate of the metals depends on their 
ionization state, which in turn is controlled by the ionizing backgrounds (UV and X-ray 
photons or cosmic rays) that were not fully considered in earlier studies. 

Here we show that the transition between the above star formation modes is driven 
primarily by fine-structure line cooling of singly-ionized carbon or neutral atomic oxygen. 
We refine earlier estimates of Z CT[t which did not explicitly distinguish between different 
coolants, by introducing separate critical abundances for carbon and oxygen, [C/H] crit and 
[0/H] cr i t , respectively, where [A/H]= log 10 (iVA/iVH) — log 10 (iVA/./VH)©, and a subscript '0' 
denotes solar values. Why are C and identified as the key species responsible for the 
global shift in the star formation mode? We find that under the temperature and density 



- 3- 



conditions that characterize population III star formation, the most important coolants are 
01 and CII whose fine-structure lines dominate over all other metal transitions 18 . Cooling 
due to molecules becomes important only at lower temperatures, and cooling due to dust 
grains only at higher densities 18,7,11 . 

The conditions that terminate the predominance of massive, population III stars, have 
dramatic implications for the reionization history of the Universe. Massive stars are much 
more efficient at producing ionizing radiation than low mass stars, due to their high surface 
temperature (~ 10 5 K). In particular, metal-free stars with a mass ^ 100M Q produce ~ 
5 x 10 4 hydrogen ionizing photons per baryon incorporated into them, a yield higher by an 
order of magnitude than obtained for a present-day mass function of stars 19,20 . Massive stars 
can therefore reionize the universe early 21 ~~ 23 at z ~ 17, as required by the recent detection 24 
of large scale polarization anisotropies in the cosmic microwave background (CMB). A strong 
UV flux just below the Lyman limit (his < 13.6 eV) is predicted 25 to be emitted by the 
same stars that are responsible for the reionization of the IGM. This soft UV radiation can 
penetrate the neutral IGM and ionize any trace amount of neutral carbon due to its low first- 
ionization potential of 11.26 eV. Since carbon is highly underabundant, the UV background 
can ionize carbon throughout the universe well before hydrogen reionization. Oxygen, on the 
other hand, will be predominantly neutral prior to reionization since its ionization potential 
is 13.62 eV. Cooling is mediated by excitations due to collisions of the respective metal with 
free electrons or hydrogen atoms. At the low fractional abundances of electrons, x e ^ 1CT 4 , 
expected in the neutral (rapidly recombining) gas at z ~ 15, collisions with hydrogen atoms 
dominate. Even if the gas had initially been ionized to a much higher degree due to the 
impact of the supernova ejecta, it would have recombined to the level stated above within 
its free-fall time [tg ~ 5 x 10 5 (n/10 4 cm -3 )^ 1 / 2 yr]. This renders our analysis independent 
of the very uncertain hydrogen-ionizing backgrounds at high redshifts (cosmic rays or soft 
X-rays) . 
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We now derive the critical C and O abundances. Our starting point is the characteristic 
state reached in primordial gas that cools only through molecular hydrogen. For the gas 
to fragment further (instead of following a slow contraction to a single massive star 1 ' 2 ), 
additional cooling due to C II or O I is required. Fragmentation requires that the radiative 
cooling rate be higher than the free-fall compressional heating rate. We find the critical metal 
abundances by equating the two rates: A-cu,oi( n 7 T) ~ 1.5nk^T/ts, where k-Q is Boltzmann's 
constant. In evaluating this threshold condition, we solve the respective rate equations for 
CII (a two- level system) and 01 (a three-level system), including all possible radiative and 
collisional transitions 18 . We assume that the emission of the fine-structure lines proceeds 
under optically thin conditions, as appropriate for the low atomic abundances considered 
here 18 . We find the critical C and O abundances by setting n ~ n c and T ~ T c , and estimate 
the uncertainties by considering a plausible range near these values (100 K ^ T ^ 200 K). 

In Fig. 1, we show the threshold C and O abundances that enable continued fragmen- 
tation for various temperature and density values. It is convenient to plot these quanti- 
ties as functions of the Jeans mass (oc T 3 / 2 n -1 / 2 ). We find critical abundance values of 
[C/H] crit ~ —3.5 ± 0.1 and [0/H] crit ~ —3.05 ± 0.2. Strictly speaking, these threshold levels 
are the required abundances in the gas phase. If a fraction of these elements were depleted 
onto dust grains, then the total C and O elemental abundances would have to be corre- 
spondingly higher. Considering the large current uncertainties in modelling the properties 
and formation channels of dust at low metal abundances 11 , we assume the limiting case of 
zero C and O depletion onto dust when comparing our theoretical predictions with observed 
stellar abundances. 

We stress that even if sufficient C or O atoms are present to further cool the gas, there 
will be a minimum attainable temperature 26 ' 6 that is set by the interaction of the atoms 
with the thermal CMB: T CMB = 2.7 K(l + z). At z ~ 15, this results in a minimum fragment 
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mass of M* ~ 2OM (n//lO 4 cm -3 ) -1 / 2 , where rif > 10 4 cm -3 is the density at which opacity- 
prevents further fragmentation 27 . It is possible that the transition from the high- mass to the 
low-mass star formation mode was modulated by the CMB temperature and was therefore 
gradual 6 , involving intermediate-mass ('population II. 5') stars 10 at intermediate redshifts. 
This transitional population can give rise to the faint SNe that have been proposed to 
explain the observed abundance patterns in metal-poor stars 15,17 . These solar- mass stars 
themselves should have formed out of metal-poor gas at relatively late times when the CMB 
temperature was sufficiently low. 

In relating our results to theoretical nucleosynthetic yields, one has to address the degree 
of mixing, or dilution, that the SN ejecta experience. SNe at high redshifts occur in dark 
matter halos with shallow potential wells, so that the resulting blast wave is predicted 28 
to escape into the expanding IGM. The dilution mass can then be estimated as M di \ ~ 
2E S n/v 2 , where E SN is the explosion kinetic energy, and v = H(z)r is the Hubble velocity 
at a distance r. Using standard cosmological parameters and ignoring radiative losses 28 , 
we get M di i ~ 1O 7 M (E S n/1O 53 ergs) 3 / 5 [(l + z)/16}~ 3 ^ 5 . In the expected mass range of 
the pair-instability supernovae 12 (PISNe) and the corresponding explosion energies E$ N — 
10 51 — 10 53 ergs, one finds dilution masses that are more than two orders of magnitude 
larger than typical values for type II SNe in nearby galaxies 14 . However, depending on 
the exact values of the ejected mass of C and O and the respective explosion energies, a 
single SN event could already lead to [C/H] and [O/H] ratios that 'overshoot' the critical 
abundances. For example, the abundances resulting from the explosion and subsequent 
mixing of a PISN at the low-mass end 12 , are [C/H]~ —2.9 and [0/H]~ —2.4, respectively, 
roughly 0.5 dex above the critical levels. Nevertheless, a single PISN at the high-mass end 
would produce final, mixed abundances still below critical. Since the dilution mass scales 
as the number of simultaneous SNe to the 3/5 power while the ejecta mass scales linearly 
with this number, multiple simultaneous SNe within the same galaxy would be less effective 
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at diluting their metal abundances into the expanding IGM. But further dilution may occur 
during the hierarchical assembly of galaxies, as infall may add metal-free gas from the IGM 
to the enriched interstellar medium. 

In Fig. 2 we compare our theoretical thresholds to the observed C and abundances 
in metal-poor dwarf 29 and giant 30 stars in the halo of our Galaxy. As can be seen, all 
data points lie above the critical O abundance but a few cases lie below the critical C 
threshold. All of these low mass stars are consistent with our model since the corresponding 
O abundances lie above the predicted threshold. The sub-critical [C/H] abundances could 
have either originated in the progenitor cloud or from the mixing of CNO-processed material 
(with carbon converted into nitrogen) into the stellar atmosphere during the red giant phase. 
To guard against this a posteriori processing of C, dwarf stars should be preferred as they 
provide a more reliable record of the primordial abundance pattern that was present at 
birth. However, since dwarfs are fainter, another strategy would be to focus on giants which 
show evidence for weak mixing (e.g. in terms of high 12 C/ 13 C isotope ratio or low 14 N 
abundance). The current data on dwarf stars (filled symbols in Fig. 2) does not yet reach C 
and O abundances that are sufficiently low to probe the theoretical cooling thresholds. Note 
also that the extremely iron-poor star HEO 107-5240 has C and O abundances that both lie 
above the respective critical levels (Bessell, M., et al., in preparation). The formation of this 
low mass star (~ O.8M ) is therefore consistent with the theoretical framework considered 
here 11 ' 17 . 

Our model suggests a powerful new probe of the first SNe to have exploded in the 
Universe. By selecting an unbiased new sample of Galactic halo stars with C or O abundances 
close to the theoretical critical values, it should be possible to sample individual SN events 
originating in the earliest epochs of cosmic star formation. The value of [C/Fe] would be 
indicative of the particular SN type. For example, the very high value of [C/Fe]~ +4 
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observed 4 in HE0107-5240, can be explained by a core collapse SN which leaves behind 
a black hole after fallback of most of the outer stellar envelope 17 . Modern surveys select 
metal-poor stars based on their weak or absent Ca II K-lines. Our model suggests that a 
better strategy to select truly 'second generation' stars is to look for stars with the lowest 
abundances of carbon and oxygen. 
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Fig. 1. — Required carbon and oxygen abundances (relative to solar values) for cooling a 
clump of metal-poor gas faster than its free-fall (compressional heating) time. The minimum 
abundances are shown as functions of the initial Jeans mass of the clump (with higher 
abundances allowing the clump to cool and further fragment). The upper panel shows the 
case of singly ionized carbon, and the lower panel refers to neutral oxygen. The curves are 
labeled by the hydrogen number density in units of cm -3 , with the solid curve marking 
the characteristic density in pre-stellar clumps of metal-poor gas, n c ~ 10 4 cm~ 3 . Heavy 
filled dots indicate particular temperature values. The critical abundance is delineated by 
the horizontal dashed lines, with a heavy line indicating the central value and the adjacent 
light lines bracketing the uncertainty within the characteristic temperature range, T c ~ 100- 
200 K, and density, ~ 10 4 cm -3 , of pre-stellar clumps. The figure indicates the minimum 
metal abundance at which fragmentation into low-mass stars starts, but it does not provide 
information about the minimum mass of the fragments at the end of the cooling process. 
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Fig. 2. — Observed abundances in low-metallicity Galactic halo stars. For both carbon 
(upper panel) and oxygen (lower panel), filled symbols correspond to samples of dwarf and 
subgiant stars 29 , and open squares to a sample of giant stars 30 . Both data sets were obtained 
at a high signal-to-noise ratio using the UVES spectrograph on the VLT. The dashed lines 
indicate the predicted level of critical carbon and oxygen abundances from Fig. 1. According 
to our theoretical model, all low-mass dwarf stars ought to lie above the critical threshold 
level of at least one of the two elements. Indeed all stars in the existing data sets satisfy 
this condition. We highlight the location of the extremely iron-poor giant star HE0107-5240 
(marked by x), whose carbon and oxygen abundances are both above the critical levels 4 . 
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